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In this paper, we study the three-dimensional (3D) iVf-flavor CP 1 model (a set of Nf CP 1 variables) 
coupled with a dynamical compact U(l) gauge field by means of Monte-Carlo simulations. This 
model is relevant to 2D s = 1/2 quantum spin models, and has a phase transition line which separates 
an ordered phase of global spin symmetry from a disordered one. From a gauge theoretical point 
of view, the ordered phase is a Higgs phase whereas the disordered phase is a confinement phase. 
We are interested in the gauge dynamics just on the critical line, in particular, whether a Coulomb- 
like deconfinement phase is realized there. This problem is quite important to clarify low-energy 
excitations in certain class of quantum spin models. If the gauge dynamics is in the deconfinement 
phase there, spinons, which transform in the fundamental representation of the SU(iVf) symmetry, 
appear as low-energy excitations. By Monte-Carlo simulations, we found that the "phase structure" 
on the criticality strongly depends on the value of Nf. For small Nf, the confinement phase is 
realized, whereas the deconfinement phase appears for sufficient large Nf > 14. This result strongly 
suggests that compact QED3 is in a deconfinement phase for sufficiently large number of flavors of 
massless fermions. 
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I. INTRODUCTION 

Some recent experiments of strongly-correlated elec- 
tron systems indicate that the usual Fermi liquid theory 
breaks down in certain cases, and the low-energy quasi- 
excitations carrying fractional/exotic quantum num- 
ber appear there. The fractional quantum Hall ef- 
fect (FQHE) 1 is a typical example, in which composite 
fermions appear as relevant excitations. Another exam- 
ple may be quantum spin models in low spatial dimen- 
sions, which have been studied quite intensively. For cer- 
tain class of s — i anti-ferromagnetic (AF) spin models 
in two dimensions, it is argued that low-energy excita- 
tions at a quantum phase transition point are spinons 2,3 . 

For studying the above "deconfined critical point" and 
the quantum phase transition itself, gauge theory is quite 
useful. Concept of confinement and deconfinement in the 
gauge theory is suitable for understanding the change of 
particle picture happening at deconfined critical points. 

In the previous paper 3 , we showed that the phase 
transition from the Neel state to the dimer state in the 
AF magnet corresponds to a Higgs (deconfinement) to 
confinement phase transition in the simple CP 1 model. 
There we were also interested in the gauge dynamics 
at the critical point. If the three-dimensional (3D) 
Coulomb-like phase is realized there as a simple loop ex- 
pansion predicts, quasi-excitations are massless spinons. 
As the low-energy excitations are magnons (spin waves) 
in the Neel state and spin-triplet excitations in the dimer 
state, the existence of spinons at the criticality indicates 
breakdown of the traditional Ginzburg-Landau (GL) the- 
ory of phase transition. This is because the GL theory 



uses an (a set of) order parameter to describe both a 
phase transition by its expectation value and low-energy 
excitations by its fluctuations in space and time. 

To study the gauge dynamics of a class of spin mod- 
els in a more general framework, we introduced the 3D 
CP 1 +U(1) lattice gauge theory in Ref. 4 . The model con- 
tains two parameters, the spin stiffness ci and the gauge 
coupling c 2 , and describes the 0(3)(CP 1 ) and 0(4) spin 
models in the specific limits, (cq — and c 2 = 00, respec- 
tively). However, from the calculation of both instanton 
density and gauge-boson mass, we concluded there that 
the confinement phase is realized on its critical line 5 . 

In this paper, we continue to study the gauge dy- 
namics of these spin systems defined in two spatial di- 
mensions at zero temperature by generalizing the above 
3D CP 1 +U(1) model to the 3D multi-flavor CP 1 +U(1) 
model. In particular, we explore the possibility of change 
of particle picture on the criticality by controlling the fla- 
vor number JVf as an adjustable parameter. 

The rest of the paper is organized as follows. In Sect. 2, 
we explain the model and its relation to the AF Heisen- 
berg model. In Sect. 3, results of Monte-Carlo simulations 
are shown. We calculated the specific heat, the gauge- 
boson mass, and the instanton density for various values 
of Nf , and found that the deconfinement phase is realized 
on the critical line for sufficiently large Nf. Section 4 is 
devoted for conclusion. 
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II. MULTI-FLAVOR CP 1 + U(l) MODEL ON THE 
3D LATTICE 

Let us first define the model on the cubic lattice, and 
explain its relation to quantum spin models. Hereafter 
we use x as the site index and [i = 1, 2, 3 as the direction 
index. On each site x, we put iVf-flavor CP 1 variables 
z", where a is the flavor index and takes a = 1, • • • , Nf. 
z™ is a two-component complex field, 



-x2 



z xli Z x2 fc 



(2.1) 



satisfying the so-called CP 1 constraint, 

= ^ \z" a \ 2 = 1 for each x and a. (2.2) 



tt=l,2 



On each link (a;, x + fj,) we put a U(l) gauge variable, 
U XI1 = exp(i9 xfi ) [6 XI1 e (-7T, +7t)]. The action of the 
model S is given as 



(2.3) 



where ci and c 2 are real parameters of the model. It is 
obvious that the action (2.3) has a local U(l) gauge sym- 
metry as well as SU(2) and SU(iVf) global symmetries; 



U 



XfJ, 



z' xa = exp(iA x )z xa , 

U' xli = exp(iA x+ll )U xlJl exp(-iA x ), 



(2.4) 
(2.5) 



6=1 

N f 

(ZxaY = E W ^ Z L, WeSV(N { ). (2.6) 



0=1 



Hereafter we call the above SU(2)(SU(iVf )) symmetry the 
spin symmetry (flavor symmetry). The partition func- 
tion Z is given by 

Z = J [dU] u(1) [dz] CP r exp(-S). (2.7) 

There are many gauge-invariant quantities composed 
of z x and U x ^. Among them, typical local combinations 
are a set of Nf 0(3) (real three-component) spins n x , 



CK —Q. CK a CK 1 

n x ~ Z x <TZ a;! n x ' n X = ' 



(2.8) 



where a =(<J\, ct 2 , 0-3)' are the Pauli matrices. For 
Nf = 1, the model reduces in the limit of c 2 = to the 
0(3) spin model described by the field n x with nearest- 
neighbor interactions 4 . For Nf > 1, there are other local 
gauge invariant objects M x ^ , which are four- component 
0(4) vectors, 



M x ^ — ^-(z x o~\z x , z x o~2Z x , z x o~3Z x , z x z x ) , 



M 



a/3* 



M 



a/3 _ 



1, M x ° = Mf*. 



(2.9) 



M x @ are complex for a ^ (3, while M x a are real. In the 
limit of c 2 = 0, one can integrate over U xfl link by link 
to obtain 



=0 = J[dz 



Z C2 =o = I [dzjcp 1 exp 



ElogIo(7x M ) 



X . fJ, 



2 JV f 



(2.10) 



where Io(7a;^) is the modified Bessel function. We note 
that M x ° are not all independent, so one needs to include 
extra interactions associated with the change of variables 
from z xa to M x @ to treat them as independent 0(4) com- 
plex spin vectors. For finite c 2 , the model involves non- 
local and/or nonpolynomial interactions among M x @ . 

On the other hand, in the limit of c 2 — > 00, the gauge 
configuration is restricted to U x ^ — 1 up to gauge trans- 
formations (2.4). Then the model reduces to an ensemble 
of independent Nf 0(4) nonlinear sigma models, 



/ [<JR]o(4) exp(ci 



[dR] 



0(4) 



x.a k=l 



R xl + iR x2 , z x2 = R x3 + iR x4 , (2.11) 



where R x is a four-component real 0(4) vector, R x 
( T? a p a \t 

[Jx xl , • • • , n, xA j . 
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FIG. 1. Phase structure of the 3D CP 1 + U(l) model 
(Nf — 1) in the C1-C2 plane obtained by the measurement 
of the specific heat 4 . The Higgs and confinement phases cor- 
respond to the Neel and dimer state of the quantum spin 
model, respectively 



In the previous paper 4 , we investigated the phase 
structure of the Nf = 1 model and found that there exist 
two phases, the ordered phase of the symmetry (2.5) and 
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the disordered phase which are separated by the second- 
order transition line c\ = C\ c (c2)- (See Fig.l.) These 
two phases correspond to the Higgs and the confine- 
ment phases in the U(l) gauge dynamics, respectively. 
In the ordered phase c\ > c\ c , there is a nonvanish- 
ing "spin magnetization" (z x crz x } ^ 0, and as a result 
the low-energy excitations are the massless components 
of z xa , which corresponds to the spin waves in the AF 
magnets (see later discussion). On the other hand, in 
the disordered phase c\ < c\ c , the confinement phase is 
realized, and the low-energy excitations are the "spin- 
triplet" vector field which is nothing but the composite 
field, n x = z x <rz x . Just on the critical line c\ = ci c , there 
is no spontaneous symmetry breaking of the internal spin 
symmetry and z xa (a — I, 2) behave as gauge-interacting 
massless bosons. Thus one may naturally expect that a 
3D Colulomb-like phase with a potential I jr may be real- 
ized there because of the screening effect by the massless 
bosons z xa . In such a phase on the critical line, the low- 
energy excitations are to be "weakly interacting massless 
spinons" z xa . With this possibility in mind, we studied 
the gauge dynamics on the critical line, and found that 
the confinement phase is realized there. This result means 
that the CP 1 model coupled with the dynamical gauge 
field (2.3) belongs to the same universality class as the 
0(3) nonlinear a model. 

The Nf — 1 case of the CP 1 model (2.3) is known to 
be a low-energy effective field theory of the nonuniform 
s = i AF Heisenberg model on a square lattice 3,6 ' 7 , 



c - t 



(2.13) 



and the CP 1 constraint, ^Z a z xa z xa \\i^ys) = |phys), re- 
stricts the magnitude of the spin to | . It was shown that 
the above Neel-dimer phase transition is nothing but the 
transition of the CP 1 model discussed above 3 . There- 
fore our investigation on the critical behavior of the CP 1 
model 4 indicates that the quantum phase transition in 
the system (2.12) belongs to the same universality class 
with the classical phase transition in the 3D 0(3) non- 
linear a model 9 . 
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FIG. 3. Ground states and excitations of the nonuniform 
AF Heisenberg model of Eq.(2.12). Each arrow represents one 
of the two spinon states z x = (1,0)* and z x = (0, 1)'. Each 
oval represents a singlet pair of NN spins, (a) dimer state, (b) 
spin-triplet excitation in the dimer state, (c) two-spinon state 
having the energy proportional to their distance (confinement 
phenomenon), (d) Neel state with an AF long-range order. 



where j is the spatial direction index (J = 1,2), S x 
is the quantum spin operator at site x, and J X j is the 
nonuniform exchange coupling. (See Fig. 2.) The ellipses 
in Eq.(2.12) represent other multi-spin interactions. By 
varying the couplings J x j , the ground state of the Hamil- 
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FIG. 2. 2D square lattice; crosses are odd sites, and filled 
circles are even sites. Solid bonds indicate that their exchange 
couplings are stronger than those on the dotted bonds. 



(2.12) changes from the Neel state to the dimer state 8 . In 
the Neel state, there exists an AF long-range order and 
the low-energy excitations are the spin waves. On the 
other hand, the ground state of the dimer state consists 
of spin-singlet pairs on nearest-neighbor (NN) sites, and 
low-energy excitations are the spin-triplets. (See Fig. 3.) 

The quantum spin operator S x is expressed by the 
Schwinger boson, i.e., the CP 1 boson operator, as follows; 



The above result supports the traditional idea for 
quantum phase transition that a quantum system in d 
spatial dimensions belongs to the same universality class 
of a certain classical system in d + z dimensions, where z 
is called dynamical critical exponent 10 . In recent years, 
however, it has been recognized that the above idea of 
the dynamical exponent breaks down in some cases; non- 
trivial physics appears at the criticality of quantum phase 
transitions. The spinons, as it was explanied above, are 
a typical example of such interesting possiblity. 

In order to see the above interesting phenomenon of 
quantum phase transition, we shall extend the model. 
The CP^ model in the 3D continuum space-time is cer- 
tainly such a model, which can be studied by the 1/N 
expansion 11 . In the leading order of the 1/N, a nontriv- 
ial infrared fixed point appears. On this fixed point, i.e., 
on the critical point, a nonlocal term for the gauge field 
Afj,(x) like, 



(2/), 



/ 1 ■ V 

d«A v - d v A 



(2.14) 



appears in the effective action due to the vacuum polar- 
ization of the massless z x . At long distances, the above 
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term dominates the usual Maxwell term which may exist 
in the original action. From (2.14), it is straightforward 
to calculate the potential energy V(r) between the two 
charges separated by a distance r as V(r) cx 1/r. Then it 
is quite interesting to study the CP^ models on the lat- 
tice for various values of N, in particular, to investigate 
the change in their critical behavors. 

Below we shall study the multi-flavor CP 1 model (2.3) 
numerically instead of the CP^ model. The reason to 
choose the multi-flavor CP 1 model is simply a matter of 
simplicity and shorter computing time in Monte-Carlo 
simulations. In the large- N limit, it is expected that the 
both models exhibit similar behavior. 12 



III. NUMERICAL RESULTS 

In this section, we present the results of our numerical 
study of the model on the 3D cubic lattice of the system 
size N = L 3 , L = 8, 12, 16 with the periodic boundary 
condition for the flavor number Nf — 1,2,3,4,5,10,14, 
and 18. We measured the internal energy, the specific 
heat, the mass of the gauge boson, and the instanton 
density. We observed no hysteresis in the internal energy 
(S)/N. 



A. Specific heat 

We first show the results of the specific heat C = 
((S—(S)) 2 ) /N measured in order to determine the phase 
structure. 
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FIG. 4. Specific heat C vs. a at c 2 = 2.0 for the N{ = 2 
case. The system size is 8 3 ,12 3 and 16 3 . C shows typical 
behavor of the second-order phase transition. 



In Fig. 4 and Fig. 5 we present C at the gauge coupling 
c-2 = 2.0 for N{ = 2 and 14, respectively. These results 
show a typical behavor of the second-order phase transi- 
tion. To confirm the transition is of second-order, we fit 
these data by the finite-size scaling hypothesis (FSSH) 13 . 
To this end, we introduce a parameter e = (ci — cioo)/cioo 
where cioo is the critical coupling in the infinite system 
(L — ► oo). Then we assume that the correlation length 
at L — ► co scales as £ cx t~ v with a critical exponent v. 
We also assume that the maximum of C at L — > oo, Coo 



diverges as Coo oc e~ a with another exponent a. Then 
FSSH predicts that the specific heat Cl (e) for the system 
size L scales as 



C L {e) = L^^L^e), 



(3.1) 



where <j>(x) is the scaling function 13 . The scaling function 
obtained from the data in Fig. 5 is shown in Fig. 6. The 
parameters are estimated as v = 1.0, cioo = 1.01 and 
(7 = 0.20. The function 4>{x) is well determined, and 
it is obvious that the FSSH is satisfied quite well. We 
investigated the phase structure of the CP^ (N — 2, 3, 4) 
models and also Nf — 2,3 cases of the multi-flavor CP 1 
models by calculating the specific heat. We conclude that 
the phase structures of both the multi-flavor CP 1 +U(1) 
model and the CP JV +U(1) model are similar to that of 
the CP 1 +U(1) model shown in Fig.l. 
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FIG. 5. Specific heat C vs. c\ at C2 
case. 
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FIG. 6. Finite-size scaling function 4>{x) of Eq.(3.1) deter- 
mined by using C of Fig. 5. All the data of L = 8, 12 and 16 
are fitted well by the single function 4>(x). 



B. Mass of gauge boson 

Now let us turn to the gauge-boson mass. We calculate 
the gauge-invariant gauge-boson mass Mq as follows 4 ' 14 . 
To define Mq, we first introduce a gauge-invariant oper- 
ator O(x), 



4 



0(x) = X! e ^ Im {U XV U X+V ^U X+I1>V U XIJ ) 

= ^ e^v sin(— 6 XV — 9 X+[J ^ + 9 X+I1 ^ U + B XI1 ), (3.2) 



/ (. . u 



where e^ u is the antisymmetric tensor. Then we intor- 
duce the Fourier transformed field 0{xz) as follows, 



6(x 3 ) = °( X Y P1X1 



We define the gauge correlation function, 

^W = ^E(oNO(x 3 + t) 

In the continuum, Dc(t) is expected to behave as 



(3.3) 



(3.4) 



D G {t) = [ dpa- 
J P 



oc e 



o 2 + M G 



(3.5) 



Typical behavor of the correlator D@(t) is shown in Fig. 7. 
We determine Mq by fitting the data Dc{t) by the ex- 
ponential form (3.5). For practical calculations, we set 
Pi = P2 = 2ir/L. 
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FIG. 7. Gauge correlation function Da{t) of Eq.(3.4) for 
N t = 18 at ci = 0.90 and c 2 = 2.0 with L = 16. 
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FIG. 8. Gauge-boson mass Ma vs. ci at C2 = 2.0 for 
Nf = !,••• ,5. Mg do not vanish for these Nt's. 



In Fig.8 and Fig. 9, we plot M G for Nf = 1 ~ 5 and 
iVf = 10, 14, 18, respectively. From the results in Fig.8, it 
is obvious that the gauge-boson mass has the minimum 
in the region close to the phase transition point. The 
minimum of the value of Mq decreases as Nf increases 
as expected, but it is still nonvanishing. In the previ- 
ous paper , we observed similar behavior of Mq in the 
CP Ar +U(l) model for N = 1,2,3,4. 
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FIG. 9. Gauge-boson mass Ma vs. ci at C2 = 2.0 for 
Nf — 10, 14 and 18. The data show that Ma vanishes 
in the critical region for Nt = 14 and 18. The case of 
Ma < implies that the square of exponent of Eq.(3.5), 
7 = Mq + p\ + P2 is smaller than 2(2n/L) 2 and we defined 
M G = -[2(2^) 2 -7] 1/2 - 



On the other hand, M G for Nf = 10, 14, 18 in Fig.9 
shows that Mq vanishes at the criticality for Nf > 14. 
This indicates that a deconfincnemt phase is realized on 
the critical line for large Nf. Appearance of the decon- 
finement phase stems from the shielding effect by the 
massless bosons z". On the critical line, low-energy ex- 
citations are massless z x and massless gauge boson xfl . 
Furthermore, we expect that topological nontrivial exci- 
tations, i.e., instantons, become irrelevant on the critical 
line due to a large number of the massless z x . (See later 
discussion.) 

From the data of Fig.8 and Fig.9, one can locate the 
minimum value of Mq along the line C2 = 2.0 for each 
Nf. These minima are presented in Fig. 10. The mini- 
mum value of Mq seems to decrease continuously as Nf 
increases. By making the linear extrapolation of the data 
for Nf = 1 ~ 5 and 10, we estimate that Mq starts to 
vanish at Nf ~ 13.5. This value is regarded as the crit- 
ical flavor number at which the phase on the criticality 
changes continuously from the confinement phase to the 
deconfinement Coulomb-like phase. 
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FIG. 10. The minimum value of gauge-boson mass Mg for 
c-2 — 2.0 vs. the flavor number Nf. The straight line interpo- 
lates the data for Nf = 1 ~ 5, 10, which intercepts Mg = 
at Nf = 13.54. 



C. Instantons 

Instantons play an important role in compact U(l) 
gauge theories 15 . Their proliferation (condensation) 
enhances fluctuations of U(l) gauge field and induces 
the confinement phase of the gauge dynamics. In the 
present 3D case, the instantons are just the magnetic 
monopoles and their condensation puts the system into 
the "dual" superconducting phase. The dual Meissner 
effect squeezes electric fluxes one-dimensionally, and as a 
result a pair of oppositely charged point particles sepa- 
rated by a distance r have the energy propotional to r, 
i.e., they are confined. 

In order to measure the instanton density, let us define 
instanton charge as in Ref. 16 . The magnetic flux 
penetrating plaquette x + /x, x + \i + z/, x + v) is given 
as, 

(-4rr < Q x ,p» < 4tt). (3.6) 

We decompose Q XiliV into its integer part 2nn Xil _ lv 
{n x ,nv is an integer) and the remaining part Q x .^ u = 
Qx,n» (mod2?r), 

@x l( ui/ = 1irn Xtftv + Q x ,nv, (— 7r < Q x ,nv < t). (3.7) 

Physically, n Xtllv describes the Dirac string. The in- 
stanton charge Q x at the cube around the dual site 
x = x + (1 + 2 + 3)/2 is defined as 

Qx ^ ^ ] £pi>pij l 'X+p J ,vp ri X U p) 

where e^p is the complete antisymmetric tensor. Then 
it is obvious that Q x measures the total flux emanating 
from the monopolc (instanton) sitting at x. The instanton 
density p is defined as 



P = Y J \Qx\/N. (3.9) 

X 

In Fig. 11, we show the instanton density p at c 2 = 2.0 
as a function c\ for the Nf — 1 and Nf = 18 cases. As 
snapshots of instanton configurations in Ref. 4 show, some 
of instantons form pairs with anti-instantons located at 
NN sites, i.e., instanton-anti-instanton dipoles. These 
dipoles are not effective for disordering the gauge-field 
dynamics and do not contribute to confinement. In fact, 
the confinement phase of the gauge dynamics is nothing 
but the plasma phase of the instantons as first shown 
by Polyakov 17 . On the other hand, the insulating phase 
of the instantons, in which almost all instantons form 
dipoles, is the deconfinement phase of the gauge dynam- 
ics. Then the density of isolated (single) instantons is a 
physical quantity which monitors whether the system is 
in the (de)confinement phase. Therefore in Fig.ll, we 
also show the density of isolated instantons p\ s , 

p is = p - 2p dp , (3.10) 

where pd p is the density of NN instanton-anti-instanton 
dipoles defined similarly as in Eq.(3.9) (The factor 2 in 
front of pd P in Eq.(3.10) comes from the fact that a dipolc 
is composed of an instanton and an anti- instanton) . 

In Fig.ll, both p and p ls almost vanish for c\ > c\ c 
in the Nf — 18 case in which Mq vanishes at the critical 
point. However, in the Nf — 1 case, there remains a finite 
instanton density at the critical point. This result and 
the calculation of Mq indicate that the Coulomb phase 
is realized on the critical line for TYf > 14. 
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FIG. 11. Density of instantons p and that of isolated in- 
stantons p is at c 2 = 2.0 vs. ci for Nf = 1 and Nf = 18. It 
is obvious that both p and p- ls for Nf = 18 tend to vanish for 
ci > cic, whereas p and pi s for Nf = 1 remain finite at the 
critical point. 
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FIG. 12. Instanton densities p and pi s at C2 = 2.0 and 
ci = 1.0 as a function of Nf. 
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FIG. 13. Log-plot of instanton densities p and pi a at 
C2 = 2.0 and c\ = 1.0 as a function of Nf. It is obvious 
that \np, \npis oc (—Nf). 



It is interesting to see how the instanton density 
changes as a function of the flavor number Nf. Re- 
sults are shown in Figs. 12 and 13 for c\ = 0.1 and 
c-2 = 2.0. From the results, it is obvious that p scales 
as p oc e~ ANt for Nf > 5, where A is a certain constant. 
This means that the main contribution in the effective 
gauge theory comes from the vacuum polarization of z" a 
(a = 1, • • • , Nf, a = 1, 2). In fact the effective gauge the- 
ory S e s(U) obtained by integrating out the CP 1 variables 
in Eq.(2.7) has the following form, 

Z = J [dU] v{1) exp(-S eff ), (3.11) 

where 



S c fi = NfS z --| (UxuUx+v^Ux+^vUx/j, + H.c.^ , 
cxp(-5 2 ) = J [dz] exp (Zx+n,aUxA + H.c.) 



(3.12) 



that S Z (U) dominates over the single-plaquette term for 
sufficiently large Nf and it determines the constant A in 
the fitting p as actually observed in Fig. 13. 



IV. CONCLUSION 

In this paper we studied the multi-flavor CP 1 model 
in three dimensions by Monte-Carlo simulations. In par- 
ticular, we are interested in the gauge dynamics on the 
critical line which separates the Higgs (Neel) and con- 
finement (dimer) phases. On the critical line, "spinons'' 
z" are massless. Their fluctuations shield the confining 
gauge force at least partly. If the number of these spinons 
is sufficiently large, the confining forces may be com- 
pletely shielded by them and the deconfining (Coulomb- 
like) force may appear instead. By calculating the gauge- 
boson mass and the instanton density, we found that the 
Coulomb-like deconfinement phase is actually realized for 
Nf > 14. The low-energy excitations on the critical line 
are the massless "spinons" z" and massless gauge boson. 
Similar deconfinement phase is expected to appear on the 
critical points of the large- TV solution of the CP^ model. 

As far as the shielding phenomenon is concerned, 
massless fermions give a similar effect as massless 
bosons. Thus the present result indicates that the parity- 
preserving QED3 with massless four-component-spinor 
fermions should have a deconfinement phase for suffi- 
ciently large flavor number of fermions, as long as the 
chiral symmetry is not spontaneously broken to avoid 
the generation of the dynamical mass. For example, in 
perturbation theory, gauge-interacting fermions generate 
the nonlocal terms like Eq.(2.14). Recently, 3D U(l) 
gauge theories coupled with gapless matter fields have 
been studied quite intensively, in particluar, to answer 
the question whether a confinement-deconfincment phase 
transition takes place 18 . The results of the present paper 
are in agreement with those obtained in these works. 

One may wonder how the results in this paper are ap- 
plied to the dynamics of realistic quantum spin mod- 
els. The corresponding quantum model for the iVf-flavor 
CP 1 model is the SU(2)x SU(iV f ) "spin" AF magnets, 
whereas for the CP N model that is SU(A^ + 1) "spin" AF 
magnets. Unfortunately, as far as we know, there are no 
materials which have the above internal quantum degrees 
of freedom. However, as we explained in the introduction 
and also in the above, study of the strongly-correlated 
electron systems often reduces to the study of certain 
gauge models of gapless matter fields. The results in the 
present paper should give an important insight into the 
phase structures of these gauge systems and the related 
strongly-correlated electron systems. 



In Eq.(3.12), J[dz] denotes the intergal over single CP 1 
field. From the above form of S e ff(U), it is expected 
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